Demonstrating anyonic fractional statistics with a six-qubit quantum simulator 
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Anyons are exotic quasiparticles living in two dimensions that do not fit into the usual categories of 
fermions and bosons, but obey a new form of fractional statistics. Following a recent proposal [Phys. 
Rev. Lett. 98, 150404 (2007)], we present an experimental demonstration of the fractional statistics 
of anyons in the Kitaev spin lattice model using a photonic quantum simulator. We dynamically 
create the ground state and excited states (which are six-qubit graph states) of the Kitaev model 
Hamiltonian, and implement the anyonic braiding and fusion operations by single-qubit rotations. A 
phase shift of n related to the anyon braiding is observed, confirming the prediction of the fractional 
statistics of Abelian 1/2-anyons. 



Quantum statistics classifies fundamental particles in 
three dimensions as bosons and fermions. Interestingly, 
in two dimensions the laws of physics permit the existence 
of exotic quasiparticles — anyons — which obey a new sta- 
tistical behavior, called fractional or braiding statistics 
That is, upon exchange of two such particles, the sys- 
tem wave function will acquire a statistical phase which 
can take any value — hence this name. Anyons have been 
predicted to live as excitations in fractional quantum Hall 
(FQH) systems 0, 0, 1^ . Alternatively, quantum states 
with anyonic excitations can be artificially designed in 
spin model systems that possess highly nontrivial ground 
states with topological order. A prominent example is the 
Kitaev spin lattice model 0, @| , which opened the avenue 
of fault-tolerant topological quantum computing 0, S] . 

It is an important goal to manipulate the anyons and 
demonstrate their exotic statistics. Towards this goal, 
a number of theoretical schemes have been proposed, 
both in the FQH regime and in the Kitaev models 
loL 11, 12, l3, 3|- However, it has proved extremely dif- 
ficult to experimentally detect the fractional statistics as- 
sociated with anyon braiding. While recent experiments 
in the FQH systems have indeed revealed some signa- 
tures of anyonic statistics [isl . [l6| . resolving individual 
anyons remains elusive Here we take a different ap- 
proach to this challenge — exploiting the spin models to 
study the anyonic statistics. Following a recent proposal 
[13], we demonstrate the fractional statistics of anyons 
by simulation of the Kitaev model on a six-photon graph 
state. The method is to dynamically create the ground 
state and excited state of the anyonic model Hamilto- 
nian, and implement the braiding and fusion operations 
by single-qubit rotations. 

How can the statistical nature of elementary particles 
be experimental observable? First let us recall the con- 
cept of quantum statistics. The wave function of a two- 
particle system ip{ri,r2) will acquire a statistical phase 



9 upon an adiabatic exchange of two particles, that is, 
'0(r2,ri) = e*^?/'(ri, r2), where = for bosons, 6 = tt 
for fermions, and 9 can be any value (0 ^ ^ vr) for 
anyons. It can be seen that a full circulation of a particle 
around the other one is equivalent to two successive par- 
ticle exchanges [§]. After such a circulation, both bosons 
and fermions show a trivial phase {(j) = 29 = 0, 27r), but 
anyons will get an observable non-trivial phase </>. To re- 
alize this idea, we need a system where anyons can be 
created and braiding operations can be carried out ex- 
perimentally. The Kitaev model is well suited for this. 

The first Kitaev model was designed on a spin lattice 
with qubits living on the edges (see Fig. la). For each 
vertex v and face /, we consider operators of the form 
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where X (Z) denotes the standard Pauli matrix i'^z)- 
These operators Ay, Bf are put together to make up the 
model Hamiltonian 



(2) 



The ground state ]'0g) of this Hamiltonian ^ is given 
by Ay\^g) — \^g) 3.11(1 B f\ip g) = \ip g) lov &\\ vcrtlccs 
and faces. Violations of these conditions cost energy and 
generate excited states J'f/'e)- A quasiparticle is created 
on the vertex Vi (face fi) if Ay^ {Bf.) acting on the ex- 
cited state IV'e), yields an eigenvalue —1 instead of +1 
for the ground state. In Ref. [1, [l3|, the quasiparticles 
on vertices are called "electric charges" (e-particles for 
short) and those on faces are called "magnetic vortices" 
(m-particles) . It is shown that these particles have un- 
usual mutual statistical properties, as we can get a phase 
flip —1 if we move one particle around the other, which 
are thus called abelian 1/2-anyons (see Fig. lb). 
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FIG. 1: The first Kitaev model 0| and anyonic braiding op- 
erations [T^]. (a). In this spin lattice, the qubits live on the 
edges and the stabilizer operators Av, Bf ([T]) represent the 
four-body interactions as illustrated, (b). Creation of quasi- 
particles and braid of an m-particle around an e-particle. A 
pair of e-particles (m-particles) are created on vertices (faces) 
by applying a Z {X) operation on the edge qubit. The quasi- 
particles can be moved horizontally and vertically by repeated 
applications oi Z or X operations. The figure shows an ex- 
ample of how an m-particle forms a closed loop around an 
e-particle through a series of moves. 
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Recently, Han, Raussendorf and Duan exploited 
the fact that the statistical properties of anyons are man- 
ifested by the underlying ground and excited states 17 1. 
So, instead of direct engineering the interactions of the 
Hamiltonian H ^ and ground-state cooling which are 
extremely demanding experimentally, an easier way is 
to dynamically create the ground state and the excita- 
tions of this model Hamiltonian, encoding the underlying 
anyonic model in a multiparticle entangled state which 
can used to simulate the dynamics of the anyonic sys- 
tem. The quasiparticles are then defined by the nega- 
tive outcome of a stabilizer element Ay or Bf. Specif- 
icaly, as illustrated in Fig. lb, with the ground state 
{ipg) prepared, one can first create a pair of e-particles 
by applying a single-qubit Z rotation. The system wave 
function will be in the excited state |'0e)- To make frac- 
tional phase experimentally detectable in a later stage, 
we apply a rotation and get a superposition state 
(l/\/2)(|^(;) + l^e))- Then we create another pair of m- 
particles and move one of them around an e-particle along 
a closed loop, and finally annihilate the m-particles. Af- 
ter doing so, it is predicted that a fractional phase tt 
will be added to |'0e), thus the superposition state will 
become (l/V2)(|V's) - IV'e)). 

As the anyons are perfectly localized quasiparticles in 
this model Hamiltonian, a small spin lattice containing 
six qubits shown in Fig. 2a allows for a proof-of-principle 
experimental demonstration 14 1 . The Hamiltonian of 
this model is Hq = -Ai - A2 - Bi - B2 - B3 - B^, 
where Ai = X1X2X3, A2 = Bi = Z1Z3Z4, 

B2 = Z2Z3Z5, B3 = ZiZe,B4,^ Z5Z6 (the subscripts of 
the Pauli matrices label the qubits). The ground state 
\tp)e of this Hamiltonian Hq is locally equivalent to a 



FIG. 2; (a). The small Kitaev spin lattice system with six 
qubits used for demonstration of anyonic braiding operations, 
(b). The six-qubit graph state which, after Hadamard (H) 
transformations on qubit 2, 3, 4, and 5, is equivalent to the 
ground state of the system in Fig. 2a. The graph state is 
associated with a graph, where each vertex denotes a qubit 
prepared in the state (l/\/2)(j0) + |1)) and each edge repre- 
sents a controlled phase gate having been applied between the 
two connected qubits p^. [19I]. The graph state is a common 
eigenstate of the stabilizer operators gi, that is, gi\'ip}6 = 
which describe the correlation in the state, and the graph 
state is the unique state fulfilling this. Here we use the same 
label as Fig. 2a-b in ref. ^4,]. (c). Experimental set-up for the 
generation of graph state and demonstration of braiding oper- 
ations. A pulsed ultraviolet laser successively passes through 
three /3-barium borate (BBO) crystals to generate three pairs 
of entangled photons 20]. The photons a, b, c, d, and / are 
combined on the three polarizing beam splitters (PBSs) step 
by step |2l[. To achieve good spatial and temporal overlap, all 
photons are spectrally filtered (AAfwhw = 3.2nm) and de- 
tected by fiber-coupled single-photon detectors (Di, • • • , Dg). 
The detector labels correspond to the qubit labels in the text 
and in Fig. 2a-b. The coincidence events are registered by 
a programmable multichannel coincidence unit. For single- 
qubit rotations and polarization analysis, quarter wave plates 
(QWPs), half wave plates (HWPs), together with polarizers 
or PBSs are used. 



six-qubit graph state [18|, [19| , which can be represented 
by the graph as depicted in Fig. 2b. This equivalence 
follows from the fact that the operators Ai, • • • , ^4 in the 
Kitaev model can be uniquely derived from the stabilizer 
operators gi (see Fig. 2b) of the graph state. 

Now we proceed with the experiment in three steps: 
(1) create and analyze of the ground state \tp)6, (2) ver- 
ify the anyonic excitations, (3) implement the braiding 
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operations and detect the anyonic phase. We use sin- 
gle photons as a real physical system to simulate the 
creation and control the anyons. The quantum states 
are encoded in the polarization of the photons which 
are robust to decoherence. The experimental set-up 
is illustrated in Fig. 3c. We start from three pairs 
of entangled photons produced by spontaneous para- 
metric down-conversion (SPDC) [20|. The photons in 
spatial modes a-b and e-f are prepared in the states 
10+),, = {l/V2){\H),\H)j + \V)^\V)j), while those in 
mode c-d are disentangled using polarizers and then pre- 
pared in the states = {l/\/2){\H)i + \V)i), where 
H{V) denotes horizontal (vertical) polarization, and i 
and j label the spatial modes. We then pass the photons 
through a linear optics network (see Fig. 3c). A coinci- 
dence detection of all six outputs corresponds exactly to 
the ground state 



|V')6 = ^m,\Hh\Hh\H),\H),\H), 

+ \V)i\V)2\V)3\H)^\H)5\H)e 
+ \H)i\H)2\V)3\V)i\V)5\V)e 
+ \V)i\V)2\H)3\V)i\V)5\V)e). 
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To verify that the ground state \ip)6 has been obtained, 
first we experimentally measure the expectation values 
of its stabilizer operators Ai, • • • , B4. These stabilizer 
operators describe the intrinsic correlations in the state 
\ip)e and uniquely define it, thus their expectation val- 
ues could serve as a good experimental signature. For an 
ideal state |V')6, all expectation values should give -1-1. 
In our experiment however, the ground state was created 
imperfectly. Figure 3a shows the measurement results, 
with all expectation values being positive in a rang be- 
tween 0.51±0.04 and 0.74±0.03, in qualitative agreement 
with the theoretical prediction. For a more complete 
and quantitative analysis, we aim to estimate the fidelity 
of the produced state, that is, its overlap with the de- 
sired one. This quantity is given by = e{4'\pexp\ip)6, 
which is equal to one for an ideal state, and 1/64 for a 
completely mixed six-qubit state. To do so, we consider 
a special observable, which allows for lower bounds on 
the fidelity, while being easily measurable with few cor- 
relation measurements (22l . [23j . By making these mea- 
surements, we estimate the fidelity of the created ground 
state to be F^^ ^ 0.532 ± 0.041 The imperfection 

of this state is mainly caused by the high-order photon 
emissions during the SPDC and the partial distinguisha- 
bility of independent photons [ist - 

We now move to the step (2). With the ground state 
\ip)e created, we apply a Z (X) rotation on qubit 3 (4), 
creating an excited state \'ipem)6 on which a pair of e- 
particles live on the vertices vi and V2, and another pair 
of m-particles on faces /i and /a (see Fig. 2a). The Z 
and X rotations are experimentally realized using HWPs 
oriented at 0° and 45°, respectively. As discussed he- 
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FIG. 3: The measured expectation values of the operators 
Ai, ■ ■ ■ , B4, oi the ground state \'4>)(i (a) and the excited state 
ii/'em)6 (b). The excited state |'!/'em)6 has a pair of e-particles 
on 111 , 112 and a pair of m-particles on /i , /s , thus the values for 
A\, A2, B\, Bz become negative. Each expectation value is 
derived from a complete set of 64 six-fold coincidence events 
in 15h in measurement basis Z®^ or X®®. The error bars 
represent one standard deviation, deduced from propagated 
Poissonian statistics of the raw detection events. 



fore, the anyonic excitations are signaled by violations 
of stabilizer conditions, that is, Ay^\%ljem)& = — iV'ewOei 
Bf- |V'em)6 = — |V'em)6 [5| ■ Thus in our case, theoretically, 
the expectation values of Ai and A2 should become — 1 
because of the e-particles, and the same for Bi and 
due to the m-particles. To verify this, we measure the 
expectation values of the operators Ai, - ■ ■ ,^4. The re- 
sults are shown in Fig. 3b, where the values of Ai, A2, 
Bi and B^ flip compared to those shown in Fig. 3a which 
supports the presence of anyonic excitations p, [l4] . 

Now we proceed to the step (3). On the ground state 
|'0)6j first we apply a ^/Z operation using a QWP ori- 
ented at 0° on the qubit 3 of the ground state yield- 
ing a superposition state |-0s)6 = (l/v^)(l'/')6 + |'0e)6), 
where |V'e)6 is the excited state with a pair of e-particles 
on vi and V2 ■ With an X rotation on the qubit 4 we fur- 
ther create a pair of m-particles on /i and /s. Then we 
perform four X operations on the qubits 6-5-3-4 to imple- 
ment the braiding operation, that is, the m-particle on /a 
is moved around the e-particle on V2 along an anticlock- 
wise closed loop. We note that the crossing at the qubit 
3, which is unavoidable in two dimensions, is relevant for 
the unusual statistics. Finally, the pair of m-particles is 
annihilated with an X operation on qubit 4 (fusion). 

After these, if there is a fractional phase (f) acquired, the 
state |V's)6 wiU become |V'/)6 = (l/V2)(|^/')6 + e*'^|^/'e)6). 
To determine this (j), we look at the correlation mea- 
surement outcomes of the six photons where the pho- 
tons 1 and 2 are fixed at |-|-) polarization, 4, 5 and 
6 at \H) and the photon 3 is measured in the basis 
(1+) + e'"|— )) with a varying in 7r/4 steps. In this set- 
ting, the six-fold coincidence counts should follow the 
relation (7(0, a) oc l-t-sin(0— a) for the state |V'/)6j thus 
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FIG. 4: (a). Measured fringes for the state \ips}6 and \ipf)G- 
The measurements in basis (|+) + e'"|— )) are done using a 
combination of HWPs, QWP and PBS. (b). The expectation 
values of the operators Ai, • • • , _B4 of the state |V'/)6 after the 
\^Zs transformation. 



an unknown phase </), if occurs, can be revealed. Fig- 
ure 4a shows the measurement results for both the state 
\^ps)6 and |V'/)6, before and after the process of m-particle 
creation, braiding and fusion. These two curves clearly 
exhibit a phase difference of tt, confirming the prediction 
of the fractional statistics. 

For a more complete proof, we implement a y/Z trans- 
formation on the qubit 3 of the remaining state |'0/)6- 
The state \tpf)6 will be converted to \tp)e if there is a frac- 
tional phase TT, otherwise it will go to |V'e)6- To test this, 
again we measure the expectation values of the opera- 
tors Ai,- ■ ■ ,B4 of the state after the ^/Z^ transformation. 
The experimental results are shown in Fig. 4b, which are 
in agreement with that the final state is and thus 
prove the fractional phase change of = tt. Here we note 
that the facts that the present setup use free-flying non- 
interacting photons and that the timescale of the braiding 
operations is extremely small (~ picoseconds for photons 
passing through the HWPs and QWPs) implies that this 
acquired phase cannot be a dynamical phase. Moreover, 
the creating, braiding and annihilating of the m-particles 
which corresponds to the operation X^XqX^X^ do not 
give rise to a phase either, as the ground state is an 
eigenstate of this observable X4XQX5X3. Similar argu- 
ments also apply to the e-particles' case. Consequently, 
this excludes possible geometrical phases [2^ and proves 
the observed phase is purely statistical. 

In summary, we have demonstrated the creation and 
manipulation of anyons in the Kitaev spin lattice model, 
and observed the fractional statistics of the abelian 1/2- 
anyons. This has been done without generating the four- 
body interactions in the model Hamiltonian but alterna- 
tively, in an easier way — by encoding the underlying any- 
onic system on a six-photon graph state, or equivalently, 
realizing the six-qubit circuit shown in Fig. 2c of ref. 
14| . It should be noted that the absence of the Hamil- 
tonian does not prevent us from studying the topological 
and statistical properties of the anyons here; however, 



topological quantum computing in a fault-tolerant man- 
ner would eventually require such a Hamiltonian. From 
a quantum-information prospective, our experiment may 
be seen as using quantum computers which have already 
well understood physics as a tool to simulate other diffi- 
cult quantum systems. Such quantum simulators can in 
principle provide exponential speedup in the simulation 
of quantum physics 27| , and may offer a more control- 



lable and clean access to study strongly correlated behav- 
iors than natural complex solid-state systems. Possible 
future work along this line may include investigation of 
the robustness of anyonic braiding Q and other possi- 
ble phases [iil using hyper- entangled graph states [j^ . 
and demonstration of some basic elements of cluster-state 
topological quantum computing [H] . 
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